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Transcription is a major obstacle for replication fork (RF) progression and a cause of genome instability. Part of this 
instability is mediated by cotranscriptional R loops, which are believed to increase by suboptimal assembly of the 
nascent messenger ribonucleoprotein particle (mRNP). However, no clear evidence exists that heterogeneous 
nuclear RNPs (hnRNPs), the basic mRNP components, prevent R-loop stabilization. Here we show that yeast 
Npl3, the most abundant RNA-binding hnRNP, prevents R-loop-mediated genome instability. npI3A cells show 
transcription-dependent and R-loop-dependent hyperrecombination and genome-wide replication obstacles as 
determined by accumulation of the Rrm3 helicase. Such obstacles preferentially occur at long and highly 
expressed genes, to which Npl3 is preferentially bound in wild-type cells, and are reduced by RNase HI 
overexpression. The resulting replication stress confers hypersensitivity to double-strand break-inducing agents. 
Therefore, our work demonstrates that mRNP factors are critical for genome integrity and opens the option of 
using them as therapeutic targets in anti-cancer treatment. 
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Eukaryotic RNA polymerase II (RNAPIl)- driven transcrip- 
tion is a tightly regulated process that occurs coordinately 
with the processing and maturation of pre-mRNAs into 
mRNAs. The formation of an export-competent messen- 
ger ribonucleoprotein particle (mRNP) and its export to 
the cytoplasm through the nuclear pore complex (NPC) 
are also coupled to transcription (Kohler and Hurt 2007). 
A number of proteins have been described as participating 
in the integration of all these processes. This includes 
THO, a conserved protein complex containing Tho2, Hprl, 
Mftl, Thp2, and Texl in Saccharomyces cerevisiae (Chavez 
et al. 2000; Pena et al. 2012), plus additional subunits in 
Diosophila and humans (Rehwinkel et al. 2004; Masuda 
et al. 2005). Together with the RNA-binding proteins Sub2/ 
UAP56 and Yral/REFl/Aly, THO forms a larger complex 
termed TREX (Strasser et al. 2002). Null mutations of THO/ 
TREX components confer similar phenotypes of transcrip- 
tion impairment, mRNA export defects, and transcription- 
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associated hyperrecombination (Luna et al. 2008), a phenom- 
enon conserved in both Caenoihabditis elegans and human 
cells depleted of THO (Dominguez-Sanchez et al. 2011; 
Castellano-Pozo et al. 2012). THO/TREX represents a par- 
adigmatic example of the connection betvi'een RNA me- 
tabolism and genome integrity that has been extended to 
other eukaryotic factors, although the molecular basis for 
this interconnection is not completely rmderstood. 

A key feature of THO/TREX mutants is their high 
levels of transcription-associated recombination (TAR). 
This has been shown to be dependent on the nascent 
RNA and the accumulation of cotranscriptional R loops 
(Huertas and Aguilera 2003), structures formed by a 
DNA-RNA hybrid and the displaced nontemplate DNA 
strand. Importantly, the connection between transcrip- 
tion, mRNA processing, and genome instability has been 
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reported by proteins involved in different nuclear mRNA 
processes. Yeast mutants of the nuclear exosome [irpS] 
and 3' end processing factors {mal4 and inalS] are also 
affected in recombination and transcript accumulation, 
although to a lesser extent than THO mutants (Luna et al. 
2005). Furthermore, the inactivation of the ASF/SF2 sphcing 
factor in chicken DT40 and HeLa cells results in genome 
instability mediated by the formation of R loops (Li and 
Manley 2005 ). Finally, a number of RNA processing factors 
have been identified in global searches performed in yeast 
and human cells as proteins preventing different forms of 
R-loop-mediated genome instability (Paulsen et al. 2009; 
Wahba et al. 2011; Stirling et al. 2012). 

In eukaryotes, the most abundant classes of RNA-binding 
proteins are the heterogeneous nuclear RNPs (hnRNPs) 
involved in multiple steps of mRNA processing and 
export (Dreyfuss et al. 2002) and the mammalian serine- 
arginine-rich (SR) proteins (SR family), which function in 
mRNA splicing, export, quality control, and translation 
(Shepard and Hertel 2009). A common feature for both 
protein families is the presence of one or more N-terminal 
RNA recognition motifs (RRMs), while their C terminus 
is composed of arginine and glycine-rich (RGG) regions in 
hnRNPs and multiple arginine and serine-rich (RS) do- 
mains in SR proteins. The most abundant hnRNP in yeast 
is Npl3, which is also considered an SR-like protein that 
shares structural homologies with both classes of RNA- 
binding proteins. Like mammalian SR proteins, Npl3 is 
a multifunctional protein involved in many processes of 
gene expression from transcription to splicing, mRNA 
export, and translation (Bucheli and Buratowski 2005; 
Dermody et al. 2008; Kress et al. 2008). It is an important 
mediator of RNA export that shuttles between the nucleus 
and the cytoplasm (Lee et al. 1996), regulated by phos- 
phorylation via Skyl kinase in the cytoplasm and de- 
phosphorylation by Glc7 phosphatase in the nucleus, 
which is also important for interaction of Npl3 with the 
export receptor Mex67 (Gilbert and Guthrie 2004). Fur- 
thermore, Npl3 stimulates RNAPII elongation and is 
an antagonist of transcription termination (Bucheli and 
Buratowski 2005; Dermody et al. 2008). Indeed, Npl3 was 
found in a complex with RNAPII and is recruited to genes 
in a transcription-dependent manner (Lei et al. 2001). In- 
terestingly, despite the general model that RNA-binding 
proteins control mRNP biogenesis and the ability of the 
nascent RNA to form R loops, no clear data reveal a role 
for Npl3 in genome instability that supports the idea of 
a failure in proper mRNA packaging that can cause RNA- 
dependent genome instability. 

Genome instability associated with most mutations in 
RNA processing factors is manifested mainly by accu- 
mulation of double-strand breaks (DSBs), hyperrecombi- 
nation, and high levels of plasmid or chromosome loss, 
implying an increase or defect in the processing of DNA 
breaks during rephcation (Aguilera and Garcia-Muse 2013). 
Consistently, some of these mutants show different de- 
grees of sensitivity to DSB-inducing agents, such as methyl 
methane-sulfonate (MMS) or hydroxyurea (HU) (Gaillard 
et al. 2009). One particularly interesting DSB-inducing 
agent is trabectedin (ET-743), which has been shown to 



form covalent adducts with guanines located in the DNA 
minor groove in vitro (Pommier et al. 1996). The action 
mechanism of this drug is particularly intriguing, and 
different reports suggest that it might lead to the gener- 
ation of replication- and transcription-dependent DSBs 
(Herrero et al. 2006; Guirouilh-Barbat et al. 2008), imply- 
ing that it may be useful in identifying new factors 
involved in the maintenance of genome integrity during 
transcription. This is of particular interest given the po- 
tential anti-tumoral properties of trabectedin that might 
serve to define new targets in cancer treatment. 

Here, after identifying Npl3 as a key factor for trabec- 
tedin resistance in S. ceievisiae, we show that npl3A cells 
have transcription-dependent genome instability pheno- 
types that are linked to a defect in replication progression. 
Such instability is at least in part mediated by R loops. 
Genome-wide chromatin immunoprecipitation (ChlP)- 
chip analyses show that Npl3 binds to highly RNAPII- 
transcribed genes and that its absence leads to R-loop- 
dependent genome-wide replication impairment, as 
detected by the accumulation of the Rrm3 helicase at 
transcribed genes. Our demonstration that an important 
and highly abundant RNA-binding hnRNP involved in 
mRNP processing and export, Npl3, has a key role in 
preventing cotranscriptional genome instability indicates 
a coordinated connection between RNA metabolism and 
DNA damage avoidance that may provide new perspec- 
tives in anti-tumor treatment. 

Results 

Hypeisensitivity to DSB-inducing genotoxic agents 
of npl3A strains 

Using a genome-wide screen in the yeast S. ceievisiae, we 
first identified the npl3A mutant as hypersensitive to 
trabectedin (ET-743), a DNA-binding drug that has been 
proposed as causing replication- and transcription-depen- 
dent DSBs (Fig. lA; Herrero et al. 2006; Guirouilh-Barbat 
et al. 2008). Next, we examined the sensitivity of npl3A 
cells to other genotoxic agents that cause DNA lesions 
that can result in stalled replication forks (RFs) and/or 
DSBs. We found that npl3A cells were sensitive to MMS, 
phleomycin (Phi), UV radiation, and HU (Fig. lA), sug- 
gesting a role of Npl3 in avoiding the incidence of DNA 
damage, in the activation of DNA damage checkpoints, 
and/or in DNA repair. 

DSB repair in npl3A cells 

To test whether hypersensitivity to genotoxic agents of 
npl3A was due to a defect in DSB repair, we generated 
double mutants lacking, simultaneously, NPL3 and the 
key genes of the two major DSB repair pathways RAD 52 
and YKU80, which affect homologous recombination (HR) 
and nonhomologous end joining (NHEJ), respectively. As 
can be seen in Figure IB, sensitivity of double mutants 
npl3A rad52A and npl3A ykuSOA to MMS was increased 
in a synergistic manner with respect to single mutants, 
suggesting that DSBs accumulate in npl3A cells and re- 
quire HR and NHEJ for their repair. 
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Figure 1. DNA repair in npl3A cells. {A) Sensitivity of BY4741 (WT) and SC086 {npl3A] strains to trabectedin (ET-743), MMS, Phi, UV 
light, and HU. [B] Sensitivity of BY4741 (WT), SC086 {npl3A], SC063 [md52A], SC087 [md52A npl3A], SC064 [ykuSOA], and SCOBS 
[ykuSOA npl3A] strains to MMS. (C) Plasmid relegation assay. [Top] Scheme of pBTM116 and enzymes used in the assay. [Bottom] 
Percentage of plasmid repair in SC071 (WT), SC073 [ykuSOA], and SC089 [npl3A] strains. (*) P< 0.05; (***) P< 0.001 (Student's t-test). 
[D] Repair of replication-borne HO-induced DSBs by SCR. [Top] Scheme of pRS316TINV plasmid carrying two inverted leu2 repeats. 
Fragments generated by HO cleavage and Xhol Spel digestion, as detected by the LEU2 probe (represented as a black line), are indicated 
with their corresponding sizes. [Left] Physical analysis of HO-induced DSB formation and its repair kinetics in WS (WT) and WSNPL3- 
01 {npl3A] strains. [Right] Quantification of DSBs (1.4-kb plus 2.4-kb bands) and SCR (4.7-kb band) related to total plasmid DNA. 
Average and standard deviation of three independent experiments are shown. 



To assay whether Npl3 was involved in one of the two 
pathways of DSB repair, NHEJ was first analyzed genet- 
ically as the ability of npl3A cells to repair a cut plasmid. 
In this assay, the 2[jl -based plasmid pBTMl 16 was linear- 
ized by restriction enzyme digestions, generating 5' or 
3' overhang ends (EcoRI and PstI, respectively), and in- 
troduced by transformation into cells (Fig. IC). Since 
plasmid maintenance requires circularization by NHEf, 



plasmid repair efficiency was calculated as the number of 
transformant colonies obtained with linearized plasmids 
relative to the tincut circular plasmid. As expected, ykuSOA 
cells showed a strong defect in DSB repair compared with 
wild type, while npl3A cells displayed a slight defect in 
the repair of EcoRI-induced DSBs and a nonsignificant 
effect in the repair of Pstl-induced DSBs (41%-53%; P < 
0.05 and P > 0.05, respectively) (Fig. IC), indicating that 
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Npl3 does not seem to have a relevant role in NHEf that 
could explain the sensitivity of npl3A cells to genotoxic 
agents. 

In S. cerevisiae, repair of DSBs occurs predominantly by 
HR using the sister chromatid as template. To test whether 
the increase in sensitivity to DSB-inducing agents of npl3A 
cells was due to an HR defect, we analyzed the involve- 
ment of Npl3 in sister chromatid recombination (SCR), 
since SCR is the preferential and major HR repair path- 
way (Kadyk and Hartwell 1992). We used the pTINV 
plasmid-borne recombinational system previously de- 
scribed (Gonzalez-Barrera et al. 2003). In this system, 
overexpression of the HO endonuclease from the GALl 
promoter generates a ssDNA break at the 24-base-pair 
(bp) mini-HO site located in one of two leu2 inverted 
repeats of the plasmid. During replication, the ssDNA 
break is converted to a DSB (Cortes-Ledesma and AguUera 
2006), which can be repaired by SCR, generating an in- 
termediate molecule that can be visualized by restriction 
enzyme digestion followed by Southern (Fig. ID). 

We determined the kinetics of the appearance of 
replication-dependent DSBs as well as the kinetics of its 
repair at different time points after HO induction (Fig. ID). 
Quantification of the different fragments generated showed 
a similar kinetics of DSB appearance, measured as the 
percentage of the 2.4-kb and 1.4-kb bands relative to total 
plasmid DNA in wild-type and npl3A cells. The value was 
maximum at 3 h after HO induction (—30% of molecules 
were cut) and decreased to <10% later on. The kinetics 
of SCR, measured as the percentage of the 4.7-kb band 
relative to the total, was also similar in both wild-type 
and npl3A cells, reaching maximum levels (—6%) at the 
9-h time point (Fig. ID). Therefore, Npl3 is not involved 
in the repair of replication-borne DSBs by SCR. It does 
not seem to play a relevant or direct role in DSB repair, 
suggesting that the sensitivity of DSB-inducing genotoxic 
agents of npl3A cells may be a consequence of repHcation 
stress induced by lack of Npl3. 

Efficient replication progression through damaged 
DNA requires Npl3 

In addition to the repair of DSBs, HR is also important for 
the repair and restart of stalled and blocked RFs (Aguilera 
and Garcia-Muse 2013). As DSB-inducing agents such as 
MMS can generate RF stalling, we reasoned that Npl3 
could play a role in the stabilization of stalled RFs or in its 
progression through damaged DNA as a way to explain 
MMS or trabectedin sensitivity of npl3A cells. To test this 
possibility, we first monitored S-phase progression in vnld- 
type and npl3A cells in the presence or absence of MMS 
by fluorescence-activated cell sorting analyses (FACS) 
(Fig. 2A). Cells were synchronized in Gl with a-factor 
and then released into fresh medium either with or without 
0.033% MMS, a concentration that has been described to 
cause extensive fork stalling in wild-type cells (Tercero 
and Diffley 2001). We found that npl3A cells exhibited 
a delay in S-phase progression in the absence of MMS, and 
this delay was much more evident in the presence of the 
drug (Fig. 2A). 



We further analyzed completion of replication after 
treatment with MMS in wild-type and npl3A cells by 
pulsed field gel electrophoresis (PFGE). This technique 
allows the distinction of fully replicated linear chromo- 
somes, which enter the gel, from actively replicating 
chromosomes, which are unable to enter the gel due to 
the branched replication intermediates and therefore 
remain in the well. Cells were synchronized in Gl with 
a-factor, treated with 0.033% MMS for I h, and then 
released into fresh medium without MMS. In both wild- 
type and npl3A cells, DNA was separated as individual 
chromosomes prior to MMS treatment, whereas most of 
the DNA was retained in the loading well 50 min after 
MMS washing and release into S phase. Interestingly, 
whereas chromosomal DNA re-entered the gel ~2 h after 
MMS exposure in wild-type cells, most of the DNA from 
npl3A cells still remained in the well after 150 min (Fig. 
2B). This result suggests that progression of the RF is 
affected in nplSA cells, possibly due to the accumula- 
tion of replication obstacles. 

A defect in replication completion after MMS treat- 
ment could be the consequence of checkpoint activation 
in response to stalled or collapsed RFs. To test this possi- 
bihty, we wondered whether the Rad53 DNA damage 
checkpoint effector kinase was activated after MMS treat- 
ment for a longer time in npl3A as compared with wild- 
type cells. We treated wild-type and npl3A cells as before 
and visualized both phosphorylated (active) and nonphos- 
phorylated (inactive) forms of Rad53 by Western. As can 
be seen in Figure 2C, Rad53 remained active much longer 
in npl3A than in wild-type cells after MMS exposure, 
indicating that the absence of Npl3 activates the S-phase 
checkpoint in response to damage. However, after repli- 
cation stress induced with HU, which depletes cells of 
dNTPs, Rad53 activation disappeared in npl3A cells as 
quickly as in wild-type cells (Fig. 2D), suggesting that 
DNA damage but not replication stress alone is toxic in 
the absence of Npl3. 

DNA breaks and transcription-associated genetic 
instabihty in npl3A cells 

Next, we tested whether the sensitivity of npl3A cells 
to MMS and their rephcation difficulties were due to a 
higher incidence of DSBs that cotdd overcome the overall 
repair capacity of the cells. To assess this, we examined 
the formation of Rad52-YFP foci, which have been shown 
previously to be a reliable method in identifying recom- 
binational repair centers (Lisby et al. 2001 ). We found that 
the percentage of cells exhibiting Rad52 foci was signif- 
icantly higher in izpBzl than in wild-type cells (8.4% vs. 3%; 
P < 0.05) (Fig. 3A), confirming that npl3A cells accumu- 
late high levels of recombinogenic DNA breaks. 

Since uniepaired DNA breaks can result in an increase 
in chromosome loss, we analyzed stability of a centro- 
meric plasmid, pRS316, in wild- type and npl3A cells after 
growing in nonselective rich medium for the same nimiber 
of generations. Figure 3B shows a significant decrease of 
more than fivefold in plasmid maintenance in npl3A cells 
(17% vs. 94% in wild type; P < 0.00 1 ), suggesting that Npl3 
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Figure 2. Replication impairment and persistent checkpoint activation in response to MMS in npl3A cells. [A] S-phase progression of 
BY4741 (WT) and SC086 [nplSA] cells synchronized in Gl with a-factor and released in the absence [left] or presence {right] of 0.033% 
MMS. [B] Analysis of completion of replication by PFGE in BY4741 (WT) and SC086 {npl3A] strains at different time points after Gl 
synchronization, treatment •with 0.033% MMS, and release into fresh medium. (C) Western against the Rad53 checkpoint kinase in 
BY4741 (WT) and SC086 {npl3A] strains. Cells were treated with MMS as in B. (D) Western against Rad53 in which Gl-arrested cells 
were treated with 200 mM HU for 3 h, washed, and released into fresh medium. A diagram of each experiment is shown at the top of 
each panel. 



plays a general role in the maintenance of genome 
stability. 

Next, we took advantage of the chromosomal leu2-k : : 
ADE2-URA3 : :leu2-k system (Fig. 3C) to assay recombi- 
nation in npl3A cells. This system is located in chromo- 
some III and contains two 2.16-kb long direct repeats of 
the leu2-k allele separated by the ADE2 and URA3 genes. 
Deletions caused by recombination between the t'wo 
repeats give rise to the loss of the URA3 marker, which 
can be scored in 5-FOA-containing medium. Consistently, 
quantification assays of Ura~ cells revealed that recom- 
bination levels in npl3A were significantly higher than in 



wild-type cells (7.3-fold; P < 0.05) (Fig. 3C), confirming the 
important role of Npl3 in the maintenance of genome 
integrity. 

Npl3 has been described as playing an important role 
in RNAPII transcription and mRNA export to the cyto- 
plasm. Consequently, we determined whether hyperre- 
combination was linked to transcription in npl3A cells. 
For this, we used different plasmid-borne recombination 
systems based on truncated repeats of 0.6 kb of the LEU2 
gene located in direct orientation and separated by in- 
tervening sequences of increasing lengths (Fig. 3D,E). 
Recombination frequencies were measured as the fre- 
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Figure 3. Genetic instability and TAR in 
npl3A cells. [A] Spontaneous Rad52-YFP foci 
formation in W303-1A (WT) and WNPL3-1D 
{npl3A ] strains and representative microscope 
images. Average and standard deviation of three 
to four independent experiments are shown. 
[B] Percentage of centromeric plasmid pRS316 
stability in W303-1A (WT) and WNPL3-1D 
[nplSA] strains. [C] Recombination analysis 
of AYW3-1B (WT) and AYNPL3-1D [npl3A] 
strains carrying the chromosomal direct 
repeat system leu2-k::ADE2-URA3::leu2-k. 
[D] Recombination analysis of W303-1 A (WT) 
and WNPL3-1D [nplSA] strains carrying L and 
LYANS plasmid systems. (£) Recombination 
analysis in the TL-lacZ plasmid system, 
whose transcription is regulated by the tet 
promoter, in the presence (Low transcription] 
or absence (High transcription) of 5 |jLg/mL 
doxycycline in W303-1A (WT), WNPL3-1D 
lnpl3A], and U678-1C {hpilA) strains. For 
B-E, average and standard deviation of three 
to four fluctuation tests consisting of the 
median value of six independent colonies 
for each one are shown. {*] P < 0.05; {**]?< 
0.01; (***) P < 0.001 (Student's t-test). For 
recombination assays, a scheme of the sys- 
tem is shown at the top of each panel. 



quency of Leu"^ colonies. First, we analyzed the effect of 
the intervening regions located between the repeats and 
found that the hyperrecombination phenotype of npl3A 
cells was observed only when a long and G+C-rich 
sequence had to be transcribed (19-fold increase in LYANS 
vs. L system; 3.3-fold in wild type; P < 0.001 and P < 0.05, 
respectively) (Fig. 3D). These results are consistent with 
published data for other factors of the mRNP biogenesis 
pathway (Luna et al. 2008). Next, we generated a new 
recombination system in which the lacZ gene is located 
between the leu2 repeats, and transcription is regulated 
by tet promoter [TL-lacZ system). As expected, recombi- 
nation was enhanced in the highly transcribed versus 
the low-transcribed system in wild-type cells (7.8-fold 
increase; P < 0.01), npl3A cells, and hprlA cells (used as 
a positive control) (Fig. 3E). However, recombination 
levels were much higher in the mutant strains compared 
with wild type (34-fold and 29-fold increase in mutants, 
respectively; P < 0.01 and P < 0.001). Therefore, npl3A 
cells show a TAR phenotype similar to that observed in 
other mRNP biogenesis factors working at the interface 
transcription-mRNA export. 

Suppression of npl3A genetic instability by RNase HI 
oveiexpiession 

We next assayed whether the TAR phenotype of npl3A 
cells was mediated by R loops, as is the case of known 
mRNP biogenesis mutants. For this purpose, we overex- 
pressed the RNase HI gene under the tet promoter from 
a newly constructed plasmid in cells carrying the RAD52- 
YFP fusion and measured the formation of Rad52 foci. We 



observed an increase of Rad52 foci in npl3A cells (Figs. 3A, 
4A). Importantly, Rad52 foci accumulation was partially 
but significantly suppressed by RNHl overexpression 
(1.7-fold vs. fivefold increase; P < 0.01), similar to the 
hpilA cells used as positive controls (0.7-fold vs. 1.7-fold 
increase; P < 0.01) (Fig. 4A). 

Next, we did the same approach using the TL-lacZ 
recombination system (Fig. 3E) and found that hyper- 
recombination was also partially but significantly sup- 
pressed in npl3A and hpilA cells (13-fold vs. 33-fold 
increase and 45-fold vs. 99-fold increase, respectively; 
P < 0.05 and P < 0.01) (Fig. 4B). We concluded that R loops 
have a significant impact over genetic stability in npl3A 
cells. We next wondered whether overexpression of AID, 
which stimulates recombination in R-loop-forming yeast 
strains (Gomez-Gonzalez and Aguilera 2007), also af- 
fected recombination in npl3A cells. As expected, AID 
overexpression from a plasmid-borne tet promoter pro- 
duced a significant increase of Leu"^ recombinants in the 
TL-lacZ system in all strains that was higher in npl3A and 
hprlA cells (2.9-fold, 4.1 -fold, and 4.3-fold increase, re- 
spectively; P < 0.001) (Supplemental Fig. SI). Therefore, 
our results indicate that npl3A -induced genome instabil- 
ity is mediated in large part by R loops, which contain 
a ssDNA fragment in the nontranscribed strand, a target 
of AID action. 

Suppiession of npl3A phenotypes by overexpression 
of hnRNPs 

Overexpression of several hnRNPs and RNA biogenesis 
factors has been shown to suppress the transcription and 
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Figure 4. Suppression of Rad52 foci and hyperrecombination 
by overexpression of RNHl. {A] Rad52-YFP foci formation in 
W303-1A (WT), WNPL3-1D {npl3A], and U678-1C [hpzlA] strains 
carrying pCM189 {-RNHl] or pCM189RNHl {+RNH1) in the 
absence of doxycycline. Representative microscope images are 
shown. Average and standard deviation of six independent ex- 
periments are shown normalized to wild-type levels. [B] Re- 
combination analysis in the TL-lacZ system in W303-1A (WT), 
WNPL3-1D {npl3A], and U678-1C {hpilA) strains under the 
same conditions as in A. Average and standard deviation of 
three to four fluctuation tests from sbc independent colonies for 
each one are shown normalized to wild-type levels. {*] P < 0.05; 
(**) P < 0.01 (Student's t-test). 

recombination-related phenotypes of transcription-RNA 
export mutants. Thus, overexpression of SUB2, THOl, 
and NAB2 suppresses mutations of either THO or THSC 
(Rondon et al. 2010). We wondered whether overexpres- 
sion of these genes also suppressed npl3A phenotypes. 
Using the chromosomal leu2-k ::ADE2-URA3 ::leu2-k re- 
combination system, we assayed the effect of overexpres- 
sion of SUB 2, THOl, and NAB2 carried in multicopy 
plasmids in both wild-type and npl3A cells. Overexpres- 
sion of these RNA factors did not have any effect on 
recombination levels in wild-type cells. The recombina- 
tion levels of nplSA cells transformed with the empty 
vector were higher than in wild type. Interestingly, over- 



expression of SUB2, THOl, and NAB 2 suppressed this 
hyperrecombination phenotype either completely or par- 
tially (7.3-fold, 3.6-fold, and 2.2-fold decrease, respec- 
tively; P < 0.01-0.05) (Fig. 5A). We then analyzed whether 
overexpression of these proteins also suppressed the tem- 
perature, HU, and MMS sensitivities of npl3A by serial 
dilution experiments. A complete or nearly complete 
suppression of the three phenotypes was observed by 
overexpressing the three genes SUB2, THOl, and NAB2 
(Fig. 5B). Altogether, these results are consistent with the 
idea that Npl3 is an mRNP biogenesis factor with a 
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Figure 5. Effect of multicopy SUB2, THOl, NAB2, and NPL3 
overexpression over different strains. [A] Recombination analy- 
sis of AYW3-1B (WT) and AYNPL3-1D {npl3A) strains carrying 
the chromosomal direct repeat system leu2-k::ADE2-URA3:: 
leu2-k and multicopy plasmids YEp351 (empty vector, YEp), 
YEpSUB2 iSUB2], YEpTHOl [THOl], and YEpNAB2 {NAB2]. 
Average and standard deviation of three to four fluctuation tests 
from six independent colonies for each one is shown. [B] Suppres- 
sion of temperature, HU, and MMS sensitivity phenotypes by 
YEpSUB2, YEpTHOl, and YEpNAB2 in serial dilutions of 
AYW3-1B (WT) and AYNPL3-1D {npl3A] cells. (C) Recombi- 
nation analysis of AYW3-1B (WT), AYNPL3-1D {npl3A], and 
AYW3-3C [hpilA] strains carrying the chromosomal direct repeat 
system leu2-k::ADE2-URA3::leu2-k and multicopy plasmids 
YEp351 (empty vector, YEp) or YEpNPL3 (NPL3). Average and 
standard deviation of six fluctuation tests from six independent 
colonies for each one are shown. (*) P < 0.05; (**) P < 0.01 
(Student's t-test). 
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function related to those hnRNPs involved in mRNP 
assembly and export. 

Reciprocally, we analyzed the effect of multicopy NPL3 
on recombination in the same chromosomal system. We 
found that multicopy NPL3 increased recombination in 
wild- type cells 10-fold [P < 0.01 ), whereas in npl3A cells, it 
did not enhance recombination further, suggesting that 
Npl3 levels have to be tightly regulated (Fig. 5C). Interest- 
ingly, overexpression of NPL3 in hpilA cells showed a 
significant reduction of hyperrecombination (5. 3 -fold de- 
crease; P < 0.05) (Fig. 5C), indicating that NplS can partially 
compensate the loss of Hprl, which is consistent with the 
idea that NplS is a structural component of the mRNR 

Genome-wide distiibution of Npl3 is eniiched toward 
the 3' end of tiansciibed genes 

To gain insight into the NplS function all over the 
genome, we performed ChlP-chip experiments using an 



NplS-MYC fusion protein. Data were subjected to com- 
putational analysis to obtain a genomic map distribution 
(Fig. 6A; Supplemental Fig. S2A) that was compared with 
previously published Hprl -Flag ChlP-chip data and used 
as a control of a cotranscriptional mRNP biogenesis 
factor (Gomez-Gonzalez et al. 2011). Statistical analysis 
revealed that 82.1% of NplS clusters and 8S.8% of Hprl 
peaked inside ORFs. Whereas NplS was significantly 
enriched in 4769 ORFs, Hprl was found in 4S24, most 
of them being coincident (7S%) (Fig. 6B; Supplemental 
Figs. S2A, S5), indicating that NplS binds to the same 
genes as THO does, in general, and that they correspond 
to actively transcribed genes. Additionally, we analyzed 
the structural and functional features of NplS -bound 
genes and found that they were longer and more ex- 
pressed than the genome average (Fig. 6B). Therefore, 
NplS binds to actively transcribed chromatin, consistent 
with its role in mRNA biogenesis. 
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We then analyzed the distribution of Npl3 along the 
length of all ORFs by subdividing them into 10 segments 
independent of ORF size, including additional upstream 
and downstream segments of the same length. The per- 
centage of clusters mapping on each segment along a 
given ORF determined the occupancy of Npl3. Interest- 
ingly, Npl3 binding increases gradually toward the 3 ' end 
of genes (Fig. 6C) regardless of gene length (Supplemental 
Fig. S2B) and similar to Hprl and other factors previously 
shown to be involved in transcription elongation (Kim 
et al. 2010; Mayer et al. 2010; Gomez-Gonzalez et al. 
2011). The result suggests that Npl3 binds during tran- 
scription, reaching a higher accumtdation as either the 
nascent RNA grows and/or the RNAPII approaches the 
3' end pause. 

Genome-wide DNA replication impaiiment in npl3A 
cells detected by Rrm3 distribution 

Results obtained so far suggest that the absence of Npl3 
causes a defect in mRNP biogenesis at the site of tran- 
scription that may cause DNA replication impairment 
and a DNA damage response deregulation responsible for 
the different forms of genome instability observed. To 
assay whether NPL3 depletion has any effect on replica- 
tion progression, we performed ChlP-chip with an Rrm3- 
Flag fusion protein. Rrm3 is a helicase required for the 
progression of the RF through obstacles in the DNA, and 
its accumulation at specific DNA sites has been used to 
identify RF pauses or stalls, which have been shown to 
occur at a high proportion in transcribed DNA regions in 
wild-type cells (Ivessa et al. 2003; Azvolinsky et al. 2009). 
In addition, we showed that Rrm3 is highly accumulated 
at transcribed genes in hprlA cells (Gomez-Gonzalez 
etal. 2011). 

We found that clusters of Rrm3 accumulation were 
distributed all over the genome in nplSA cells (Supple- 
mental Fig. S3) and that such clusters were significantly 
longer in nplSA cells (average size of 656 bp vs. 475 bp in 
wild type; P < 0.01), with a total extension of 3.12 Mb vs. 
2.50 Mb in wild-type cells (26% and 21% of the genome, 
respectively) (Fig. 7 A], consistent with the idea that 
replication obstacles are stronger or more abundant in 
nplSA cells. Indeed, the percentage of clusters mapping at 
ORFs was also higher in the mutant (83.6% vs. 76.0%, 
respectively), consistent with the effect on replication 
being preferentially observed at transcribed regions and 
with the role of Npl3 promoting RNAPII transcription. 
Furthermore, both clusters of Rrm3 and Npl3 were sig- 
nificantly enriched at introns, snRNAs, snoRNAs, and 
RNAPIII- transcribed genes, including tRNAs (Supplemen- 
tal Fig. S4; Supplemental Table SI), suggesting an addi- 
tional role of Npl3 facilitating RF progression through 
these regions. Additionally, we found that Rrm3 bound to 
centromeres in npl3A cells, but not Npl3. 

There was a 52% coincidence between the ORFs re- 
cruiting Rrm3 in wild-type and npl3A cells, the ORFs 
showing a high overlap with Npl3-bound ORFs (Supple- 
mental Fig. S5). To determine the main functional and 
structural features of the genes to which Rrm3 preferen- 



tially binds, we selected the top 500 ORFs with the 
highest signal log2 ratio average in a sliding window of 
200 bp. Statistical analysis revealed that, on average, 
these were long, GC-rich, and highly transcribed genes in 
both wild-type and npl3A cells, indicating that RNAPII- 
transcribed genes are more prone to affect RF progression 
(Fig. 7B). Interestingly, top Rrm3-bound genes were more 
expressed in npl3A cells, suggesting that replication ob- 
stacles occur preferentially at the highest transcribed 
genes in the absence of Npl3 (Fig. 7B), these genes pref- 
erentially being shorter all over the genome (Marin et al. 
2003). In this sense, although gene ontology analysis 
revealed that genes involved in ribosome biogenesis 
or intron-containing genes were significantly overrepre- 
sented in npl3A cells, this is consistent with the fact that 
ribosome biogenesis genes are highly transcribed and 
represent the major class of intron-containing genes in 
yeast. 

It is worth noting that Rrm3 -binding profiles showed 
an increase toward the 3' ends of genes in both wild-type 
and npl3A cells, with a high decrease in the regions 
immediately upstream of and downstream from them 
(Supplemental Fig. S6A), consistent with the pattern of 
Npl3 binding to transcribed genes (Fig. 6C). We also found 
a slight displacement to the 5' end of long genes in the 
Rrm3 profile in the mutant that correlates with the Npl3 
profile (Supplemental Figs. S2B, S6B), suggesting that 
Npl3 helps in preventing or resolving RF progression 
impairment all over the gene, especially at long genes. 
Additionally, we measured the Rrm3 signal log2 ratio 
average along all ORF segments and found a notable 
increase in npl3A cells (Fig. 7C), which indicates that 
replication obstacles are stronger in RNAPII genes when 
Npl3 is absent. 

Finally, microarray analysis revealed that gene expres- 
sion profiles of wild-type and npl3A cells were quite 
similar (Pearson's R correlation coefficient of 0.978), 
indicating that the absence of Npl3 does not have a 
significant impact on global expression levels. Neverthe- 
less, down-regulated genes in npl3A cells were longer, 
G-i-C-richer, and better expressed than up-regulated genes 
(Supplemental Fig. S7A), consistent with the results ob- 
tained in hpilA, tho2A, andsub2zl cells (Gomez-Gonzalez 
et al. 2011). Furthermore, we found a good correlation 
between Rrm3 binding and expression levels in the 2758 
genes accumulating Rrm3 in both the wild-type and the 
npl3A cells in the mutant (Pearson's R correlation co- 
efficient of 0.59234) that was clearly higher than that of 
wild-type cells (0.24348, Supplemental Fig. S7B), support- 
ing the idea that replication obstacles occur preferentially 
at highly transcribed genes in npl3A cells. 

Rrm3 accumulation at transcribed ORFs in npl3A cells 
is R-loop-dependent 

We showed that R-loop formation at transcribed genes is 
in part responsible for the genetic instability phenotypes 
of npl3A cells. If this instability were mediated by RF 
progression impairment, we would expect that R loops 
could also contribute to the replication impairment 
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Figure 7. Rrm3 genome-wide recruitment in wild-type and npl3A cells. (^) Genomic view of Rrm3 recruitment in WRBb-9B (WT) and 
N3RBb-4B {npl3A] cells. A fragment of chromosome II is plotted with the signal log2 ratio values. Red histograms represent the 
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across the average ORF plotted as signal log2 ratio average per each segment. (D) Detailed analysis of ChlP-chip data of Rrm3-Flag at the 
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caused by the absence of Npl3. To test this, we performed 
specific Rrm3-Flag ChIP RT-qPCR experiments in se- 
lected genes with high levels of Rrm3 recruitment as 
determined by the ChlP-chip, such as GCN4 (Fig. 7D). As 
can be seen in Figure 7E, there was a significant enhance- 
ment of Rrm3 recruitment in this gene in npl3A cells, 
validating the ChlP-chip results. Importantly, overex- 
pression of RNase HI significantly reduced the levels of 
Rrm3 to those of the wild type [P < 0.05), suggesting that 
the genome-wide effect of Npl3 in preventing transcription- 
dependent rephcation impairment is linked to cotranscrip- 
tional R-loop formation. 



Discussion 

In this study, we show that the RNA-binding hnRNP Npl3 
has a function in preventing transcription-dependent 
genome instability. Moreover, loss of Npl3 results in 
increased sensitivity to DSB-inducing agents and ge- 
nome-wide replication impairment mediated by the for- 
mation of R loops. Our results not only provide a new link 
between mRNP biogenesis and the prevention of R-loop- 
mediated transcription-replication collisions but support 
a model in which a number of RNA-binding hnRNPs 
prevent RNA-mediated replication stress and open the 
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possibility of using them as a therapeutic target in anti- 
cancer treatment. 

We identified npl3A as a mutant hypersensitive to 
trabectedin (ET-743). The action mechanism of this anti- 
tumor drug remains poorly understood, but it has been 
shown that in vitro, it forms covalent adducts with 
guanines located in the DNA minor groove (Pommier 
et al. 1996). Addition of trabectedin leads to the genera- 
tion of DSBs, the activation of the DNA damage check- 
point, and the recruitment of proteins of the nucleotide 
excision repair and HR pathv\fays in Schizosacchaiomyces 
pombe (Herrero et al. 2006; Guirouilh-Barbat et al. 2008). 
Interestingly, it is worth noting that trabectedin binding 
to the DNA leads to the formation of a complex that 
resembles a DNA-RNA hybrid (Marco et al. 2002). npl3A 
cells were also hypersensitive to other DNA-damaging 
agents, including MMS, Phi, HU, and UV light, all of 
which lead directly or indirectly to RF stalling and/or 
DSBs. However, the ability of npl3A cells to repair DSBs 
by SCR or NHEJ is not affected or is poorly affected 
(Fig. 1), suggesting that the lack of Npl3 generates replica- 
tion stress that results in the accumulation of recombino- 
genic DNA breaks, which is strengthened by genotoxic 
agents. In fact, sensitivity of nplSA cells to MMS was 
increased synergistically in npl3A rad52A and npl3A 
ykuSOA cells, in agreement with a previous genome- vdde 
analysis that reported a synthetic growth defect for npl3A 
md52A double mutants (Pan et al. 2006). Interestingly, 
npl3A cells are not able to complete replication after 
treatment with MMS, leading to activation of the DNA 
damage response to prevent the collapse of putatively 
stalled RFs. Altogether, these data suggest that npl3A 
cells are under permanent replication stress. 

In the absence of induced DNA damage, npl3A cells 
accumulate DNA breaks monitored as Rad52 foci and 
show plasmid instability and transcription-dependent 
hyperrecombination (Fig. 3), which is consistent with the 
identification of Npl3 in a screening for mutants involved 
in yeast minichromosome maintenance (Wahba et al. 
201 1). Importantly, hyperrecombination in the absence of 
Npl3 is transcription-dependent, according to the role of 
Npl3 in transcription (Lei et al. 2001 ), and can be partially 
suppressed by the overexpression of RNase HI, which 
can also suppress Rad52 foci accumulation (Fig. 4), in- 
dicating that R loops are at least in part responsible for 
the genetic instability generated in npl3A cells. Addition- 
ally, human AID overexpression exacerbates recombina- 
tion, consistent with its enhanced action on the displaced 
ssDNA of cotranscriptional R loops (Chaudhuri et al. 
2003; Gomez-Gonzalez and Aguilera 2007). R loops form 
naturally during specific cellular processes, including 
Escherichia coh plasmid replication, mitochondrial DNA 
replication, and immunoglobulin class switching (Aguilera 
and Garcia-Muse 2012). However, they can also be gener- 
ated as transcriptional by-products that compromise ge- 
nome integrity, as shown for cells (from yeast to humans) 
depleted of a number of conserved proteins involved in 
RNA metabolism, including the THO-TREX complex 
(Huertas and Aguilera 2003; Dominguez-Sanchez et al. 
201 1; Castellano-Pozo et al. 2012), the DNA-RNA heU- 



case Senl/Senataxin (Mischo et al. 2011), topoisomerase I 
(Tuduri et al. 2009; El Hage et al. 2010), or the ASF/SF2 
spHcing factor (Li and Manley 2005). In addition, several 
screenings in yeast and human cells have revealed that the 
genome instabihty generated by the absence of a great 
nimiber of RNA processing factors is also R-loop-dependent 
(Paulsen et al. 2009; Wahba et al. 2011; Stirling et al. 2012). 

An enhancement of R-loop formation or stabilization 
might explain the transcription-dependent and R-loop- 
dependent replication collisions and genome instability 
of npl3A cells, consistent with its main putative role in 
mRNP assembly, although we cannot discard an indirect 
contribution of a defect in mRNA export given the in- 
teraction of Npl3 with the Mex67 export factor, whose 
mutations also lead to genetic instability (Gilbert and 
Guthrie 2004). Importantly, however, overexpression of 
SUB2, THOl, and NAB2 is able to suppress npl3A pheno- 
types (Fig. 5), suggesting an interchangeable role of these 
proteins that rely on their ability to bind nascent RNA. 
Consistent with the functional interrelationship between 
Npl3 and THO, Npl3 overexpression is able to partially 
suppress the phenotypes of hptlA. Indeed, the fact that 
overexpression of NPL3 leads to hyperrecombination in 
wild- type cells, as has been also shown for SUB2 over- 
expression (Fan et al. 2001), suggests that the proper 
stoichiometry of proteins involved in mRNP assembly is 
key to preventing RNA processing defects from compro- 
mising genome integrity. However, we cannot discard 
that the inhibition of transcription termination reported 
by excess of Npl3 (Bucheli and Buratowski 2005) could be 
one of the causes of this instability. 

Many transcription elongation and RNA processing 
factors, including THO-TREX, have been shown to be 
recruited along transcribed ORFs in a gradient profile that 
increases toward the 3' end of genes all over the genome 
(Kim et al. 2010; Mayer et al. 2010; Gomez-Gonzalez 
et al. 2011). This, together with the 3' end processing 
defects of THO mutants (Saguez et al. 2008) as well as the 
transcription termination role of SenI RNA-DNA heli- 
case (Kawauchi et al. 2008), suggests that termination 
may be one of the most sensitive steps connected to 
RNA-mediated genome instability. Notably, Npl3 is also 
distributed along transcribed ORFs according to an in- 
creasing pattern toward the 3' end (Fig. 6) and causes a 
preferential reduction in RNA accumulation of long and 
highly transcribed genes in concordance with the features 
of the genes to which it is preferentially bound. These 
results confirm that Npl3 has a main role in cotranscrip- 
tional mRNP assembly. 

Given the double effect of npl3A on gene expression 
and genome instability, the main question remaining is 
whether both are interconnected. It is known from many 
reports from bacteria to human cells that transcription 
stimulates recombination (Gaillard et al. 2013). Consis- 
tent with recombination being the main pathway of 
repair of DSBs generated during replication, it is believed 
that most TAR events arise during replication as a result 
of transcription-replication collisions that cause the stall- 
ing and collapse of forks (Azvohnsky et al. 2009; Bermejo 
et al. 2009). Here we show that, indeed, npl3A cells have 
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a replication delay and that most replication obstacles 
occur at transcribed genes in these cells, visualized as the 
genome-wide distribution of Rrm3, a DNA helicase re- 
quired for RFs to bypass obstacles (Ivessa et al. 2003; 
Azvolinsky et al. 2009). As Npl3-bound genes, Rrm3- 
bound genes are also long and highly expressed (Fig. 7), 
and there is a high enrichment of Rrm3 in ribosomal 
protein genes, the most highly transcribed class of genes 
in yeast. This is consistent with a previous study showing 
that Npl3 is preferentially associated with ribosomal 
protein mRNAs and other highly transcribed transcripts 
(Kim Guisbert et al. 2005) and the reported role of Npl3 
mediating nuclear export of large ribosomal subunits 
(Hackmann et al. 2011). Rrm3 is preferentially observed 
in a gradient manner that increases toward the 3' end 
of RNAPII genes, following a pattern similar to that 
previously shown for Rrm3 in THO mutants (Gomez- 
Gonzalez et al. 2011), which is partially suppressed by 
RNase HI overexpression (Fig. 7E). Therefore, cotran- 
scriptional R loops in npl3A cells constitute a major cause 
of RF stalling responsible for recombination-mediated 
genome instabiUty. 

Consequently, we propose a model in which Npl3 acts 
as a key player in coupling transcription of RNAPII genes 
with mRNP formation and export, thus facilitating rep- 
hcation through highly transcribed regions (Supplemen- 
tal Fig. S8). In the absence of Npl3, the mRNP assembly 
would be suboptimal, and the nascent mRNA may hy- 
bridize more stably with the transcribed DNA strand, 
forming R loops that would constitute roadblocks to rep- 
hcation and cause fork stalling and breaks that would 
demand HR for fork restart and survival. In summary, 
our results show that Npl3, the most abundant RNA- 
binding hnRNP, is necessary to prevent R-loop-mediated 
transcription-rephcation conflicts and genome instabil- 
ity. Given the relevance of genome instability as a com- 
mon feature of tumoral cells (Bartkova et al. 2005; GorgouUs 
et al. 2005), the R-loop-like structure of the trabectedin- 
DNA complex (Marco et al. 2002), and the importance of 
Npl3 for resistance to this anti-tumor drug, it would be 
interesting to explore whether deregulation of RNA- 
binding proteins leading to the formation of R loops in 
human cells could be used to identify new anti-tumoral 
targets in cancer treatment. 

Materials and methods 

Yeast strains and plasmids 

Yeast strains and plasmids used in this study are listed in 
Supplemental Tables S2 and S3. 

Physical analysis of SCR 

SCR assays were carried out essentially as described (see the 
Supplemental Material; Gonzalez-Barrera et al. 2003). 

Chip analyses 

Recruitment of Rrm3 to chromatin was determined by CHIP 
analyses on the endogenous GCN4 gene as previously described 
(see the Supplemental Material; Hecht et al. 1999). 



Micioaiiay gene expression analysis 

Microarray determination of total RNA was performed using the 
Affymetrix platform (see the Supplemental Material) as previ- 
ously described (see the Supplemental Material; Gomez-Gonzalez 
et al. 2011). The expression data can be accessed at Gene Ex- 
pression Omnibus (GSE50187 and GSE50186). 

ChlP-chip experiments 

S. cerevisiae ohgonucleotide tiling microarrays were provided by 
Affymetrix. The high-density ohgonucleotide arrays used were 
able to analyze yeast chromosomes at a 300-bp resolution, each 
of the 300-bp regions being covered by at least 60 probes. ChlP- 
chip of asynchronously growing cells was carried out as described 
(see the Supplemental Material; Katou et al. 2006; Bermejo et al. 
2009). The ChlP-chip data can be accessed at Gene Expression 
Omnibus (GSE50187 and GSE50185). 

Statistical analysis of genome-wide data 

Microarray data were normalized by RMA (robust microarray 
average) and statistically analyzed by LIMMA (linear models for 
microarray analysis), comparing the mutant expression profile 
with its isogenic wild-type strain. The genes showing at least 
a 1.5-fold expression change with a P-value < 0.01 with a false 
discovery rate (FDR) correction were considered as altered. ChlP- 
chip data were analyzed using the Tiling Array suite (TAS) soft- 
ware from Affymetrix. For each probe position, TAS produces the 
signal and the change P-value, taking into account the probes 
localized within a given bandwidth around the inspected probe. 
Protein chromosomal distribution was then analyzed by detect- 
ing binding clusters, which were defined as ranges within the 
chromosome respecting the following conditions: estimated 
signal (IP/SUP-binding ratio) positive in the whole range, P-value 
< 0.01, minimum run of 100 bp, and maximum gap of 250 bp. 
The results were visualized with the University of California at 
Santa Cruz Genome Browser, developed and maintained by the 
Genome Bioinformatics Group (Center for Biomolecular Science 
and Engineering at the University of California at Santa Cruz; 
http://genome.ucsc.edu). For statistical analysis of the func- 
tional and structural features of the genes, expression levels 
were taken from microarray of wild-type ceUs. Distribution of 
binding sites along genes was carried out as previously described 
(Gomez-Gonzalez et al. 2011). 

Miscellanea 

Analysis of sensitivity to genotoxic agents, PFGE, Western 
analyses, Southerns, FACS using a FACScahbur Becton Dickinson 
machine, and yeast cultures were performed using standard 
procedures. Genetic assays of NHEJ, plasmid maintenance, 
recombination, and Rad52 foci are described in the Supplemen- 
tal Material. 
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